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Abstract 

Using supersymmetric quantum mechanics we develop a new me- 
thod for constructing quasi-exactly solvable (QES) potentials with 
two known eigenstates. This method is extended for constructing 
conditionally-exactly solvable potentials (CES). The considered QES 
potentials at certain vàlues of parameters become exactly solvable and 
can be treated as CES ones. 
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1 Introduction 



Since the appearance of quantum mechanics there has been continual inter- 
est in models for which the corresponding Schrodinger equation is exactly 
solvable. With regards to solvability of the Schrodinger equation there are 
three interesting classes of the potentials. 

The first class is the exactly solvable potentials allowing to obtain in ex- 
plícit form all energy levels and corresponding wave functions. The hydrogen 
atom and harmònic oscillator are the best-known examples of this type. 

The second class is the so-called quasi-exactly solvable (QES) potentials 
for which a finite number of eigenstates of corresponding Hamiltonian can 
be found exactly in explicit form. The first examples of QES potentials were 
given in [1-4]. Subsequently several methods for generating QES potentials 
were worked out and as a result many QES potentials were found [5-13] (see 
also review book [Q). Three different methods that are based respectively 
on the polynomial ansatz for wave functions, the point canonical transfor- 
mation, and the supersymmetric (SUSY) quantum mechanics are described 
I2| . Recently, an anti-isospectral transformation called also as duality 



m 



transformation was introduced in [pT5|| . This transformation relates the en- 
ergy levels and wave functions of two QES potentials. In |Tü| a new QES 



potential was discovered using this anti-isospectral transformation. 

The third class is the conditionally-exactly solvable (CES) potentials for 
which the eigenvalues problem for the corresponding Hamiltonian is exactly 
solvable only when the parameters of the potential obey certain conditions. 
Such a class of potentials was first considered in fT7fl . It is interesting to note 
that in |ï|] it was demonstrated the equivalence of the condition required for 
potential obtained in W?\ to be a CES potential with the condition that this 
potential can be put in an explicitly supersymmetric form. Recently, new 
examples of CES potentials have been discovered [TT5], |2Ü], pT| . 

Very useful algebraic tool for investigation of the problem of exact solv- 
ability of the Schrodinger equation is the SUSY quantum mechanics (for 
review of SUSY quantum mechanics see [^2], ^3|). For constructing QES po- 
tentials the SUSY method was used in [10-12]. The starting point of this 
method is some initial QES potential with n + 1 known eigenstates. Then 
applying the technique of SUSY quantum mechanics one can calculate the 
supersymmetric partner of the QES potential which is a new QES potential 
with n known eigenstates. In the papers |20, |2Ï| the SUSY quantum me- 
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chanics was used to develop some generalized method for the constructing 
the CES potentials. 

In our previous paper []23| we have proposed a new SUSY method for 
constructing QES potentials with two known eigenstates which in contrast 
to the papers [10-12] does not require knowledge of the initial QES potentials 



for a generation of a new QES ones. In |25| we extended this method for 
constructing QES potentials with three known eigenstates. 

The present paper is devoted to further development of the SUSY method 
proposed in j24| and to extension of this method for constructing CES poten- 
tials. In result we obtain new QES and CES potentials. A new interesting 
point is that QES potentials with two known eigenstates become exactly solv- 
able at certain fixed vàlues of parameter and therefore they can be treated 
as CES potentials. 



2 Supersymmetric quantum mechanics 

In the Witten's model of supersymmetric quantum mechanics the SUSY 
partner Hamiltonians H± read 
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where 



V±(x) = i (W 2 (x) ± W'(x)) , W'(x) = ^^-, (3) 

W(x) is referred to as a superpotential. In this paper we shall consider the 
systems on the full real line — oo < x < oo. 

The eigenvalues and eigenfunctions ip^{x) of the Hamiltonians H± 
are related by SUSY transformations which in the case of unbroken SUSY 
read 

K + x = K, E~ = 0, (4) 
ChO) = -T^B + ip+ (x), ip+(x) = r — B~i;- +1 (x). (5) 

E-n+l 
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As a consequence of SUSY the Hamiltonians H + and í/_ have the same 
energy spectrum except for the zero energy ground state. The latter exists 
in the case of the unbroken SUSY. Only one of the Hamiltonians H± has a 
square integrable eigenfunction corresponding to the zero energy. Here we 
use the convention that the zero energy eigenstate belongs to H_. Due to 
the factorization of the Hamiltonians H± (see ([I])) the ground state for H 
satisfies the equation B~ip^{x) = the solution of which is 

i/>o(x) = C exp (- J W(x)dx^ , (6) 

Cq is the normalization constant. Here and below C denotes the normal- 
ization constant of the corresponding wave function. From the condition of 
square integrability of the wave function i/jq(x) it follows that the superpo- 
tential must satisfy the condition 

sign(ÏU(±oo)) = ±1, (7) 

which is the condition of the existence of unbroken SUSY. 

For a detailed description of SUSY quantum mechanics and its application 
for the exact calculation of eigenstates of Hamiltonians see reviews [P2[ , f23|l . 
Just the properties of the unbroken SUSY quantum mechanics reflected in 
SUSY transformation (£|), (d) can be used for exact calculation of the energy 
spectrum and wave functions. In the present paper we use these properties 
for the generation of the QES potentials with two known eigenstates and 
CES potentials. 



3 QES potentials with two known eigenstates 

We shall solve the eigenvalue problem for the Hamiltonian H_. The ground 
state of this Hamiltonian is known and is given by wave function (|GJ) with 
energy Eq =0. In order to calculate the excited state of H_ we use the 
following well-known procedure used in SUSY quantum mechanics. Let us 
consider the SUSY partner of H_, i.e. the Hamiltonian H + . If we calculate 
the ground state of H + we immediately find the first excited state of H 
using the SUSY transformations (^), (|5|). In order to calculate the ground 
state of H + let us rewrite it in the following form 

H+ = H W + e = B+Bf + e, e > 0, (8) 



4 



which leads to the following relation between the potential energies 

V+(x) = V<: 1 \x) + e, (9) 

and superpotentials 

W 2 (x) + W'(x) = W?(x) - W[{x) + 2e, (10) 

where e is the energy of the ground state of H + since we supposed that 
similarly to H_ has zero energy ground state, Bf and Vl 1 \x) are given by 
(0) and (^) with the new superpotential W\(x). 

As we see from the ground state wave function of H + is also the ground 
state wave function of íf_ and it satisfies the equation Bï%Pq{x) = 0. The 
solution of this equation is 

xp+(x) = C+ exp (- J W!(x)dx\ , (11) 

where for the square integrability of this function the superpotential W\(x) 
must satisfy the same condition as W(x) (|7j). Using (||) and (||) we obtain the 
energy level Eï = e and the wave function of the first excited state ipï(x) 
for 

Repeating this procedure in the case of shape invariant potentials [26| and 



self-similar potentials [57], |2^| it is possible to calculate all energy spectrum 
and the corresponding wave functions. As a result for these cases many 
exactly solvable potentials were obtained |29] (see also review |f22f ). 

We consider a more general case and do not restrict ourselves to the shape 
invariant potentials or self-similar potentials. In this case it is not possible 
to obtain all energy spectrum. In |2j| we obtained the general solution of 
equation (|ÏÜ1) and thus derived a general expression for QES potential with 
two explicitly known eigenstates. The bàsic idea consists in finding such a 
pair of W(x) and W\{x) that satisfies equation (|Ï0P . To do this we rewrite 
equation (TO) in the following form 

W' + {x) = W-(x)W+(x) + 2e, (12) 

where 

W+(x) = Wi(x)+W(x), (13) 
W_(x) = W x {x)-W{x). 



5 



This new equation ([12]) can be easily solved with respect to W-(x) for a 
given arbitrary function W + (x) or with respect to W+(x) for a given arbitrary 
function W-(x). Then from ( P^ ) we obtain superpotentials W(x) and Wi(íc) 
which satisfy equation (^D[). 



In contrast to our paper [24| where we use the solution with respect to 
W-{x), in the present paper we use the solution with respect to W+(x). 
This solution is more convenient for construction CES potent ials and can be 
written in the following form 



W+(x) = exp (^J dxW-(x)\ 2e J dxexp (— J dxW. 



x) + A 



(14) 



here A is the constant of integration. 



In order to simplify solution (JJ) let us choose W-(x) to be of the form 

W-{x) = -cj)"(x)/(j) , (x). (15) 

To provide a nonsingularity of W-(x) and as a result nonsingularity of V±(x) 
we shall consider a nonsingular monotonic function 4>(x) satisfying the con- 
dition 4>'{x) > 0. Then substituting (|Ï5|) into (|ï^ ) we obtain 



= (2e0(x) + \)/4>\x). 



(16) 



Note that the constant À can be included into the function <p(x) and thus 
for W+(se) we obtain 

W+{x) = 2t(j>{x)/(j>'{x). (17) 
Finally, for superpotentials W(x) and W\{x) we have 



W(x) = (e<f>(x) + -cf>"(x))/<P'{x), 
Wx(x) = (e<f>{x) - l -(j>"{x))/ct>'{x). 



(18) 
(19) 



Using this result for the wave functions of the ground state with the 
energy Eq = and excited state with Eï = e we obtain 



-1/2 



exp 



dx(j){x) / '</>' (x) 



Eo = 0, 



-1/2 



exp 



dx(j){x)/(t)'{x) , E{ 



€. 



(20) 
(21) 
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Note that as we see from ( |T3| ) W+(x) must satisfy the same condition 
(0) as W(x) and W\(x) do. Then because (f>(x) is monotonic function and 
4>'(x) > from ( |TT| ) it follows that has one node. Therefore ipï{x) given 
by (|2Ï|) also has one node and thus corresponds to the first excited state. 
The functions <fi(x) that satisfy described condition provide also the square 



integrability of the wave functions (|20| ) and (|2T|). 

It is worth to stress that <fi{x) = ipï (x) / ipQ (x) from which follows an 
interesting fact. Namely, the ratio 4>{x) of the wave functions of first excited 
state and ground state and the distance e between the corresponding energy 
levels entirely determine the potential energy. 

QES potential V_(x) is given by @ with superpotential (|T8|) . Choosing 
different <p(x) and e we obtain different QES potentials with explicitly known 
two eigenstates. 

Now let us consider interesting examples of new QES potentials which 
become exactly solvable at certain fixed vàlues of parameter e and thus can 
be treated as CES ones. 



3.1 Example 1 

Let us put 

<f ) {x) = i3H 2k+1 {ix), (22) 

where H m (x) is Hermite polynomial, the final result does not depend on the 
constant f3. The superpotentials in this case read 

W(x)= 1 x + i2k( 1 + l) H * k - l{Í *\ (23) 

H 2k {ix) 

W 1 (x) = 7 * + t2k ( 7 - 1) ^ff , (24) 
where we have introduced the notation 

^ÏTT < 25 > 

Substituting the superpotential W(x) (^3|) into (||) we obtain the following 
QES potential V-(x) 
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V-{x) = + 2k(2k - 1) ( 7 + l) 2 H * k ~f x) (26) 

2 H 2k (ix) 

-2F( 7 + l)( 7 + 3)(^iM) 2 
\ H 2k {ix) ) 

+kj (7 + 1) - ^7. 

The wave functions of the ground and first excited states read 

il,o (x) = C (iJ 2fc (zx))- (1+ ^ 2 exp(-7x 2 /2), (27) 
^-(x) = C^H 2k+l {ix) (H 2k {ix))-^ )/2 exp(- 7 a; 2 /2). (28) 

It is interesting to note that in the special case 7 = 1, i.e. 

e = 2k + l, (29) 

the second term in Wi(x) drops up and W\(x) corresponds to the superpo- 
tential of a linear harmònic oscillator. Then VÍ (x) and, as a result of @, 
V + (x) are the potential energies of the linear harmònic oscillator. Therefore 
in this case the SUSY partner H + is the Hamiltonian of the linear harmònic 
oscillator and we know all its eigenfunctions in explicit form. Using SUSY 
transformations (f|), @ we can easily calculate the energy levels and the 
wave functions of all the excited states of H . The energy spectrum of H 
in this special case is the following 

E =0, E~ = n + 2k, n = 1,2, ... (30) 

Thus QES potential (|26| ) at fixed value of e (|29| ) becomes exactly solvable 
and therefore can be treated as CES potential. Note that V-(x) in this special 
case 7 = 1 corresponds to the potential obtained by Bagrov and Samsonov 



[ p0| , |3T| via Darboux method and latter by Junker and Roy [£(], ^TJ within 
SUSY approach. 

3.2 Example 2 

Consider the function 

,/ \ H 2k+ i(ix) 

<j){x) =13 TT-r- , k>m, 31 
H 2m (ix) 
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which generalizes the one given in the first example. For superpotentials we 
obtain 

W(x) = x • • i\m-^~ 



H 2m (ix) 



_ % e + 2k - 2m + 1 

H 2m+1 (ix) / H 2m (ix) - H 2k+2 (ix)/ H 2k+ i(ix)' 



WAx) = x + iAm- 



H 2m (ix) 



i e-(2fc-2m + l) 

H 2m+ i(ix) / H 2m (ix) - H 2k+1 (ix) / H 2k (ix) 

The QES potential VI (x) is given by (§) with superpotential (0). The 
expressions for the QES potential and the wave functions in this case are 
somewhat complicated and we do not write them down in explicit form. 
We would like to stress the following very interesting point. In the special 

case 

e = 2k - 2m + 1 (34) 

the second term in superpotential W\ (x) ( |33|) drops up and then this W\ (x) 
coincides with superpotential W(x) ( p3|) from the first example for 7 = 1 
which is exactly solvable. Then using the same explanation as in the end of 
first example we may conclude that for the potential VI (x) calculated with 
superpotential (|32|) in the special case ([34]) it is possible to get all energy 
levels and corresponding wave functions and thus V_(x) can be treated as 
CES potential. For the energy levels we obtain 

E = 0, Eï = 2k-2m+l, 
E~ = n + 2k, n = 2,3,4,... (35) 

In this special case the potential energy V_ (x) coincides with the one studied 
in i|. 



4 CES potentials 

In this section we develop a subsequent method for constructing the CES 
potentials using the results of the previous section. 
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Suppose that Wi(x) is a given superpotential that corresponds to the 
exactly solvable potential V_(x). The example of such superpotential is a 
shape invariant one |p6| |. As a result of (H) V + (x) is also exactly solvable and 



thus for H + we know in explicit form all energy levels and the corresponding 
eigenfunctions. Then using SUSY transformations (Q), (|5|) we can easily 
calculate all excited energy levels and wave functions of its SUSY partner 
H^, the wave function of the ground state is given by (|6|). But to do this we 
must have the superpotential W(x) which is expressed through Because 
Wi(x) is a given function it is convenient to represent the superpotential 
W(x) using ( |Ï8|) and (|l9|) in the following form 

W(x) = W 1 (x) + (36) 

Similarly, the new exactly solvable potential V^(x) can be written as follows 

= ÏMM + WI( X) ) + (f||) 2 + wwjg - «. (37) 

In this expressions the function is not arbitrary one but must satisfy 
( [19]) for a given Thus we must solve equation (|T9| ) with respect to 

for a given Wi(sc) which can be written in the following form 

^<p"{x) + W^x^'ix) = e<f)(x). (38) 

In order to transform this equation to Schròdinger type equation let us 
write 4>(x) in the form 

0( x ) = f( x ) exp f- J dxWx(x) \ . (39) 
The new function f(x) satisfies the equation which can be rewritten as follows 
-\f"(x) + vl 1) (x)f(x) = -ef(x), (40) 

where 

VÍ 1 \x) = \{Wl{x) + W[{x)). (41) 
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As we see it is a Schrodinger type equation of SUSY quantum mechanics 
but with negative energy. The sign in the right hand side of equation (fíO|) 
can be changed using a duality transformation called also as anti-isospectral 
transformation [15 

i = ix. (42) 

Then equation fl4~0f) reads 



i d 2 no 



2 de 

where we have introduced the notations 



■V. 




/HO, 



(43) 

(44) 
(45) 
(46) 

In the present paper we shall consider only such superpotentials W\(x) for 
which the dual superpotential is real function of £. Then equation (pj|) 

is ordinary Shròdinger equation of SUSY quantum mechanics. Using (|39|) 
and ([H]) the solutions of equation ([38]) can be expressed via the solutions of 
equation (H3) in the following form 




x 



7(0 exp /de^i(e) 



f(ix) 
fo(ix) 



(47) 



where /(^) is the solution of equation ([43]) for the energy e, /o(£) is the 
solution of the same equation for the zero energy e = 0. Note, that for 
obtaining <p{x) it is not necessary to use only square integrable solutions of 
equation fl43|). The solutions must be such that 4>(x) is a monotonic function 
that has one node. 

Thus, now we have the problem to solve equation (f43|). In order to solve 
this equation the following fact is important. If W\(x) is such that the 
corresponding SUSY partners V±\x) belong to the class of shape invariant 
potentials then W\(£) gives also the shape invariant SUSY partners V± (£). 
To see this recali that the superpotential in the shape invariant case satisfies 
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the following equation [2B| 



W^x, a) + = W?(x, ai) ~, ; + 2i2, (48) 

where the superpotential W\{x) in the left and right hand sides of this equa- 
tion have different vàlues of parameters a and a%, the remainder R does not 
depend on x. From this equation using the definition ([46]) we obtain 



W? (6 ai) + ^ J = a) + 2i2. (49) 

As we see also satisfies the shape invariant equation. Note, that in 

compare to fl48|) the set of parameters a in the equation (f49[) is replaced by 
«i and vice versa. 

Thus gives the shape invariant SUSY partners V±\í) and equation 

( pO] ) can be solved exactly. Using these solutions on the basis of ( pT7| ) we get 
0(x). Substituting <f>(x) into fl3"7|) we obtain a new exactly solvable potential 
V-(x) which is lower SUSY partner of the known exactly solvable potential 
V+(x). Of course we must verify that <p(x) satisfies the considered above 
conditions, namely, <j)(x) must be monotonic function with one node. 

In conclusion of this section let us consider some examples. 



4.1 Example 3 

Let us put 

W x {x) = x. (50) 

Such a choice corresponds to linear harmònic oscillator. The dual superpo- 
tential in this case has the same form as W\{x) 

Wi£) = tï (51) 

and thus equation ( |4"3"D is the Shròdinger equation for linear harmònic oscilla- 
tor. Then using the well-known wave functions of stationary states of linear 
harmònic oscillator on the basis of ( f£7| ) we obtain <p(x) = 0H2k+i(ix) (P is 
some constant), which is exactly equal to (0). Here each k corresponds to 
the energy e = 2k + 1. Note, that in order to satisfy appropriate conditions 
for <p{x) we select only odd solutions of equation (|43|). 

This example reproduces the result of the first example of section 3 in 
the special case 7 = 1. 
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4.2 Example 4 

Let us consider a superpotential W\(x) that corresponds to the Rosen-Morse 
oscillator. In this case 

W x (x) = atanh(x) (52) 
is shape invariant. The dual superpotential 

Wi(£) = atan(0 (53) 
is also shape invariant. The potential energy corresponding to (|53|) reads 

lÇJ 2cos 2 (0 2 1 J 

As we see the dual potential energy has singularities at the points 

£ n = ! + 7m, n = 0,±l,±2,.... (55) 



Traditionally, Schròdinger equation (Ï3) with potential energy ( |51| ) is con- 
sidered on the interval between two neighbouring singularities using zero 
boundary conditions for solutions /(£)• Without any problem we consider 
the solutions of ( f43| ) on the full £-line which take zero value in all points (|55|) 

/(e«) = 0, n = 0,±l,±2,.... (56) 

Such solutions can be easily obtained with the help of SUSY quantum me- 
chanics. Using three first odd solutions of ( fj3|) for 0(x) given by ( pTFD we 
obtain 

4>i(x) = sinh(x), for e = e\ (57) 
03 0*0 = [1 — a + (2 + a) cosh(2x)] sinhx, /or e = €3, (58) 
05 (x) = [6 + a + 3a — 4(a 2 + 2a - 3) cosh(2x) + 

(a 2 + 7a + 12) cosh(4x)] sinhx, for e = e 5 , (59) 

where 

e k = ((a + k) 2 -a 2 )/2. (60) 

Note that we may directly verify that functions <f>(x) given by fl57|) , (|58|) , fl5"9|) 
indeed satisfy equation (|38|) . Substituting the obtained 4>k(x) and e = into 
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fl37|) we get the set of exactly solvable potentials V~(x, k) (here we explicitly 
write down the dependence of the potential on k) 

V_{x, k) = tanh 2 (x) (^-a(a - 1) + <f> k (x)($ k (x) + 2a)) - e k + |, (61) 



where 



$ 3 (x) 



$i(x) = 1, 
3(2 + a) cosh(2x) + a + 3 



$ 5 (x) 



(2 + a) cosh(2x) - a - 1 
(3 + a)[5(4 + a) cosh(4x) + 4(5 - a) cosh(2a;)] + a(5 - a) + 30 
(3 + a)[(4 + a) cosh(4x) - 4(1 + a) cosh(2s)] + 3(1 + a)(2 + a) ' 



The potential VL(j;, 1) reproduces the Rosen-Morse one. Other potentials 
are new exactly solvable ones. 

Note that for k = 3 the potential can be written in the following explicit 
form 

V(x o)- 4 ( 3 + 2Q ) i 4 ( 1 + Q ) (62) 

l ' ; ((2 + a)cosh(2x)-l-a) 2 (2 + a) cosh(2x) - 1 - a y ' 

(l + a)(2 + a) (3 + a) 2 



2 cosh 2 x 



and was previously obtained by us in |^4| . 

Note, that e is the parameter of the superpotential and thus the parameter 
of the potentials V±(x). As we see V-(x) is exactly solvable when e is equal 
to a certain fixed value. It is worth to stress that for given functions <f)(x) 



(for example (|58D , (p9|)) and arbitrary e it is always possible using the results 



of previous section to construct QES potentials with two known eigenstates. 



These QES potentials at certain fixed vàlues of e ( 60|) become exactly solvable 
(|6"ID and can be treated as CES potentials. 



5 Conclusions 

We have developed a new SUSY method for constructing QES potentials 
for which we know in explicit form energy levels and wave functions of the 
ground and first excited states. From the obtained general expressions for 
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QES potential and wave functions of the ground and first excited states the 
following interesting fact can be derived. The ratio of the wave functions of 
first excited state and ground state and the distance between corresponding 
energy levels entirely determine the potential energy. 

The method developed for constructing QES potentials with two known 
eigenstates is extended for generating CES potentials which are exactly solv- 
able at certain fixed vàlues of parameter e. Finally, this new exactly solvable 
potential is the lower SUSY partner to the known exactly solvable poten- 
tial. In this sense our method for generating CES potentials is similar to the 
method proposed in [[ÏÏ], ^TJ although the realization is different. In addition, 
our method gives the interesting relation between QES and CES potentials. 
Namely, when the parameter e of QES potentials is equal to a certain fixed 
vàlues then QES potentials become exactly solvable and can be treated as 
CES ones. 

Note, that important moment in our approach for generating CES po- 
tentials is the duality transformation which we use to transform equation 
( |38|) to Schròdinger type equation. In the present paper we consider only 
superpotential for which the dual one ( |46|) is a real function of new variable. 
The case of complex dual superpotential is more complicated and we plan 
to consider this case in future. It will give a possibility to extend the class 
of CES potentials which can be obtained by suggested in the present paper 
method. 
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